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The Equation of State and the properties of matter in the high temperature deconfined 
phase are analyzed by a quasiparticle approach for T > 1.2 Tc- In order to fix the parameters 
of our model we employ the lattice QCD data of energy density and pressure. First we 
consider the pure SU(3) gluon plasma and it turns out that such a system can be described in 
terms of a gluon condensate and of gluonic quasiparticles whose effective number of degrees 
of freedom and mass decrease with increasing temperature. Then we analyze QCD with 
finite quark masses. In this case the numerical lattice data for energy density and pressure 
can be fitted assuming that the system consists of a mixture of gluon quasiparticles, fermion 
quasiparticles, boson correlated pairs (corresponding to in-medium mesonic states) and gluon 
condensate. We find that the effective number of boson degrees of freedom and the in- medium 
fermion masses decrease with increasing temperature. At T ~ 1.5 Tc only the correlated pairs 
corresponding to the mesonic nonet survive and they completely disappear at T ~ 2 Tc- The 
temperature dependence of the velocity of sound of the various quasiparticles, the effects of 
the breaking of conformal invariance and the thermodynamic consistency are discussed in 
detail. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Nq 
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I. INTRODUCTION 



The state of matter at asymptotic high temperature is predicted unambiguously by Quantum 
Chromodynamics (QCD). Indeed, it is by now weh-estabhshed that, at sufficiently high tempera- 
ture T, the hadronic constituents form a deconfined and chirally symmetric Quark-Gluon Plasma 
(QGP). In this regime, due to the asymptotic freedom of QCD, the QGP is expected to be a weakly 
interacting gas of quark- and gluon-quasiparticles. 

However, the various heavy-ion experiments at Brookhaven, GSI and CERN laboratories, de- 
voted to the creation and detection of new forms of highly excited matter indicate that, at the 
reachable energy scales, the produced phase exhibits a strong collective behavior which is incom- 
patible with a weakly interacting QGP. 

In particular, standard perturbative QCD cross sections for quarks and gluons do not allow for 
a rapid thermalization and fail to reproduce, by hydrodynamical models, the observed phenomena 
(e.g. the so-called elliptic flow). Moreover, the results obtained at the Relativistic Heavy Ion 
Collider (RHIC) seem to show that, at large temperature and small baryon density, the new state 
of matter can be described as an almost perfect fluid of strongly interacting particles. 

At present the best available tool to study the non-perturbative properties of QCD in these 
extreme conditions is the numerical simulation on the lattice. Lattice QCD (IQCD) [l|, [3] gives 
precise indications on the deviations of the pressure, p, and of the energy density, e, of the system 
from the Stefan-Boltzman (SB) values of an ideal gas of quarks and gluons even at temperatures 
as high as T ^ 4 Tc- This suggests that remnants of the confining interaction play a non-trivial 
role in the deconfined phase [3|]. 

In order to motivate our approach to the deconfined phase of QCD let us first consider how 
the properties of matter in the confined phase change with increasing temperature. It is well 
known that the hadronic content of matter and its dynamical properties change dramatically with 
increasing temperature [1, d, @, Q]- In particular, for T = 100 MeV, the largest fraction of particles 
(~ 82%) are pions, a smaller fraction, ~ 12%, consists of strange mesons and only about 6% of 
particles are in excited states. The mean free path can be estimated to be ^ 15 fm and the 
volume available per particle is roughly 30 fm^. In these conditions hadrons do not overlap and the 
system can be described as a mixture of the various hadronic states. With increasing temperature 
the density increases and the hadronic content changes. As an example, for T = 150 MeV, the 
excited particles are ~ 40%, the remaining 60% of particles are pions and strange mesons; the 
mean free path is ~ 2fm and the volume per particle is 4 fm^. Therefore in the confined phase, 
with increasing temperature, the percentage of particles in excited states increases and close to the 
transition temperature it is reasonable to expect that the density of particles is large with a large 
fraction in excited states. Since the mean free path can be of the order of the average distance 
between particles the system may have liquid like properties. 

For temperature larger than the critical temperature, it is not clear which are the correct degrees 
of freedom. In any case since the lattice simulations show that the transition from the confined 
phase to the deconfined phase is a smooth crossover @, @|, we expect that for T > Tc mesonic states 
may survive in the deconfined phase. This conclusion is confirmed by the surprising results of Refs. 
T) U) l3] where, by evaluating the static two point correlation function of a quark-antiquark pair. 



it has been shown that correlated qq states survive up to temperatures ~ 2Tc. Moreover it seems 
that the dissociation temperature for the correlated states depends on the mass of the quarks. 



t—L 



the higher the mass the higher the melting temperature 13|, llJ]. This implies that a detailed 
analysis is required to understand the particle content of matter or, more generally, the effective 
degrees of freedom above Tc. The presence of mesonic states in the deconfined phase, observed in 
IQCD, has also been discussed by various authors in the random-phase approximation and by 
employing phenomenological models 16, 17, 18, [l^ and has been shown ^] that may lead to a 
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rapid thermalization of heavy quarks in the QGP. 

Since correlated mesonic states may survive in the deconfined phase they have to be included in 
a consistent thermodynamical description of the system above Tc- In our previous paper [21], we 
have, indeed, shown that IQCD data of pressure, p, and energy density, e, in the range T = 1.2 — 2 
Tc, are consistent with a description based on q,q,g quasiparticles, relatively light boson states 
and gluon condensate. 

In the present paper we shall review in more details those results and analyze in a more general 
way the thermodynamics of the system and the role of the gluon condensate. Due to the short 
mean free path, a quasi-particle approach is probably not reliable for temperature very close to the 
transition point, therefore we expect our results to be valid for temperatures larger than ~ 1.2 Tc- 

In Section II we discuss the gluon plasma i.e. pure SU(3) gauge theory. The lattice results of 
the so-called interaction measure, indicate deviations of the system from the conformal symmetric 
behavior. We find that part of this deviation is due to the gluon condensate. However, the gluon 
condensate is not enough to fit the lattice values of the trace anomaly at finite temperature, that 
measures the breaking of the conformal symmetry. Therefore there must be other mechanisms that 
break the scale invariance of the system (dynamically) generating a new scale. In Section III Al 
where we discuss a quasiparticle model for the gluon plasma, we assume that such a mechanism can 
be described in terms of massive gluon quasiparticle whose effective number of degrees of freedom 
is also temperature dependent. Since at asymptotic high temperatures QCD is approximately scale 
invariant we consistently find that the effective gluon mass decreases with increasing temperatures. 
In Section IIIII we extend our analysis to the Quark-Gluon Plasma and show that a quasiparticle 
description can reproduce the lattice results of the energy density and of the pressure of the QGP. 
This in turn leads to a constraint on the number of effective boson degrees of freedom and on the 
mass of the fermion quasiparticles. Section [IV] is devoted to conclusions and outlooks. 



II. ANALYSIS OF THE GLUON PLASMA 



Lattice simulations of pressure and energy density provide useful information regarding the 
properties of QCD in a large range of temperatures (for a review see Ref. [22.]). One interesting 
point is to understand whether the system can be approximately scale invariant. We know that 
QCD at zero temperature is not a conformal invariant theory: the coupling constant depends on a 
typical energy scale, Aqcd, and the stress energy tensor Q^^'^ has a trace anomaly. In this Section 
we analyze the pure SU(3) gluon plasma, where it is assumed that the masses of the fermions are 
infinitely large, i.e. in the so-called quenched approximation; the analysis of full QCD with finite 
quark masses (unquenched approximation) will be developed in the next Section. 

The trace of the average energy momentum tensor is related to the energy density, e, and the 
pressure, p, of the system by the relation: 

e;i = e-3p. (1) 



At finite temperature, one can write [23| 



e^(r) = e(T) - 3p{T) =<G^ >o- <G^ >T, (2) 



where = —gG'Jj^^G^'^ and < G^ >o and < G^ >t are the gluon condensate at zero and at finite 
temperature respectively. Note that the energy density e(r) and pressure p{T) are consistently 
normalized to zero at T = 0. 

Equation ([2|) states that the breaking of the conformal invariance at finite temperature is only 
due to the gluon condensate. We can check this relation by using different lattice results. Indeed 
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in Ref. [2J, |25l] the gauge invariant two-point correlation function of the gauge field strengths have 
been evaluated on the lattice by using the correlator method. In particular in Ref. f25| the gluon 
condensate has been evaluated at finite temperature in both pure gauge and full-QCD. It turns 
out that, at moderate temperatures T ^ 2Tc, the chromo-magnetic component < of the 

gluon condensate survives above the deconfining temperature and is (within the statistical errors) 
temperature independent: 



< G2 >^ 



1, (3) 



whereas the ratio between the chromo-electric contributions < G^ at zero and at finite temper- 
ature, 

< G2 >g ^ ' 

rapidly decreases above the deconfining temperature. 

Since a T = the electric and magnetic terms are equal, defining, 

Ai(T) =-ij<G2>o[l-cr(r)], (5) 

and assuming that in a pure gluon plasma the conformal invariance is only broken by the gluon 
condensate, one has 

^ = Ai(T). (6) 

Now considering the lattice results of energy density, e^, and pressure, pL, of Ref. HI, 3 we can 
obtain the "interaction measure" for a pure SU(3) gluon plasma, 

A(r) = iiffl_|5i(I), (7) 

and we can compare such expression with Ai(T). 

In Fig. [T]we show the plots of A(T), full line (black online), and of the Ai(T), dashed line (red 
online). In this plot we are assuming that the gluon condensate at zero temperature is given by, 
< >o~ 0.03 GeV^ 0, li^, Q and c§"fr} is obtained by interpolation of the data obtained in 



lattice QCD simulations reported in Ref. [25 1 



It is clear from Fig. [T] that the gluon condensate is not able to describe the trace anomaly and 
other contributions must be present that stem from interactions and/or masses. 

The gluon condensate plays an important dynamical role in QCD. Indeed, at zero temperature 
it contributes to a negative non-perturbative vacuum energy. In general, a bosonic condensate is 
a macroscopically populated state with zero momentum and the contribution of the modes with 
/c 7^ is small. For these reasons, we will assume that at any temperature the gluon condensate 
does not significantly contribute to the pressure but does contribute to the energy density of the 
system. 

In order to consider the effect of the gluon condensate on the energy density we consider the 
quantity — Ai. In Fig. [2] we show the plots corresponding to the lattice data of pressure, pL, 
versus energy density, e^, dashed line (red online) and the data of pressure, pL, versus ei — Ai. 



1 1.5 2 2.5 3 3.5 4 

T/T^. 



FIG. 1: (color online) Interaction measure A(r) defined in Eq.([7|), full line (black online), and contributions 
of the gluon condensate Ai, dashed line (red online), defined in Eq. ((S]) as a function of T/Tc in the pure 
SU(3) Yang-Mills. 




FIG. 2: (color online) Plots of pressure as a function of energy density. The diagonal line (green online) 
corresponds to the conformal solution with p = e/3. The dashed line (red online) corresponds to the lattice 
data of energy and pressure, whereas the dot-dashed (back online) line corresponds to the quasiparticle 
gluons, with pressure equal to pl and energy — Ai. The (blue) dot on the diagonal line corresponds to 
the ideal SU(3) gas solution. 

The effect of subtracting the gluon condensate is of making the relation between energy density 
and pressure of the gluon system more similar to the conformal solution, that corresponds to the 
diagonal full line p = e/3. However, subtracting the gluon condensate is not enough for making 
the system conformal invariant. Note that for p = and e = the system seems to be conformal, 
but this happens because the pressure and energy density have been normalized to be zero in the 
confined phase. The (blue) dot on the diagonal line corresponds to the energy density and pressure 
of an ideal gluon gas. 

Therefore one can conclude that the trace anomaly at finite temperature is not saturated by 
gluon condensation and that the description of the system requires other dynamical ingredients. 
The rigorous analysis in Ref.fi^, where a gluon gas with physical state space reduced to the 
so-called Gribov fundamental modular region has been considered, gives essentially the same indi- 
cation. 

In the next sections we shall consider a more phenomenological approach, based on the intro- 
duction of massive gluon quasi-particle, to describe the difference between trace anomaly and gluon 
condensate shown in Figs. [T]and[2j 
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A. Quasiparticle model for the Gluon Plasma 

Let us recall some thermodynamical relations. The free energy density of a system is given by 

f = -LlnZ{T,V), (8) 

where Z(T, V) is the partition function, and energy density and pressure are given by, 

r2 51nZ(r, V) 



V dT 



(9) 



. ^ t'-^ . (10) 



For homogeneous systems, in the thermodynamic limit, 

lnZ{T, V) dlnZ{T, V) 



(11) 



V dV 

and therefore the pressure can directly be obtained from the free energy density, 

P = -f. (12) 

Using this relation one can express the entropy density s and the trace anomaly in terms of 
derivatives of the pressure with respect to temperature, 

e + P dp 

' = -r = dT^ (13) 

<e;i>T=6-3p = T5^(p/T4). (14) 

Note that for an ideal gas p ^ T'^ and the trace of the stress energy tensor is equal to zero. 
However, when interactions and/or masses are present the pressure is not proportional to and 
one has to include these effects in the trace anomaly. 

An effective method to take into account the interaction is by considering temperature depen- 
dent gluon masses and temperature dependent gluonic degrees of freedom. Therefore in our analysis 
we work in a quasiparticle picture and all the in-medium effects are treated, at the mean field level, 
as an effective mass, Mg(T), and as an effective number of degrees of freedom, Dg{T), for the 
gluons. The values of these parameters will depend on the temperature and this dependence will 
be determined by fitting the lattice data of pressure and energy density where the gluon condensate 
contribution will be properly taken into account. However, temperature dependent parameters and 
also the presence of the gluon condensate require that the thermodynamical consistency must be 
carefully checked. This point is discussed in details in the Appendix A. 

According to the previous discussion, let us rewrite Eq. ^ for a gluon plasma as follows: 

< >T= e-3p =<G^ >o- < >T +Aeg , (15) 

where A0g represents the contribution to the trace anomaly due to gluon quasiparticles. 

As discussed, the gluon condensate is an important dynamical ingredient and we will assume 
that it contributes to the energy density but not to the pressure of the system. Therefore in our 
quasiparticle picture, where the contributions to the thermodynamical quantities come from gluon 
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condensation and gluons with in medium mass Mg(T), the pressure and energy density of the 
system are given by, 

P = Pg e = £9 + econ > (16) 

where, 

= T'A^{T) = 1 < g2 >o [1 - , (17) 

is contribution to the energy density due to the gluon condensate. The quasiparticle contributions 
to pressure and energy density are respectively, 

__log(l-e-«/^)-i = D,{T)p,{T), 

^9iT) = D,{T) /(^^)3 ,.JT_ 1 = D,{T)e,{T) , (18) 

with DgiT) the temperature dependent number of gluon degrees of freedom, Pg{T) and eg(T) the 
gluonic pressure and energy density per degree of freedom, and 

u;g{k,T) = ^k^ + MgiT)\ (19) 

is the gluon dispersion law. 



B. Determination of Mg and Dg 



In order to obtain the functions Dg (T) and Mg (T) we will employ the data of the pressure and 
energy density obtained in lattice simulations reported in Refs. [11,3. We obtain Dg{T) and Mg{T) 
substituting 

Pl{T)=p{T) eL{T) = e{T), (20) 



in the Equations ()16p and (jlSp . and solving for the effective mass of the gluons and the effective 
number of the degrees of freedom. In particular the gluonic mass can be obtained solving the 
equation, 

gL(r)-6g,",(r) _ igjT) 

Pl{T) Pg{T)' ^ ^ 

and the effective number of gluonic degrees of freedom can be obtained from the relation, 

D^^T) = llfl , (22) 

where the solution for the effective gluon mass obtained solving Eq. (pT]) has been plugged in pg (T) . 

Before presenting our results for Dg{T) and Mg{T), it is useful to clarify some points regarding 
the range of values of the temperature where our results are reliable. Let us consider the equation 

A(r) = Ai(T) + D,iT){ig^-3p,) ^ (23) 
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which at the critical temperature imphes 

D,iT^)Ce -3p,)\t. ^ ^^^^^ _ ^^^^^^ ^ ^2^^ 

From this equation and from Eq. (j2ip one can see that close to Tc the values of Mg{T) and 
Dg{T) will strongly depend on the values of the lattice data. 

Indeed, for T close to Tc there are two solutions of Eq. (I2ip : Mg(Tc) — > (corresponding to the 
saturation of the trace anomaly with the gluon condensate) but also Mg(Tc) — > oo. Moreover, for 
T — > Tc, we have that dc ^ and the actual value for Dg{T) will depend on the ratio between two 
small quantities. In particular we find that when the gluon mass diverges also Dg(Tc) diverges, 
unless dc is exactly zero at Tc. In any case, as previously discussed, a quasi-particle approach is 
unreliable at the critical point and for these reasons our numerical analysis will be limited to the 
range ~ 1.2 — 4 Tc- 

From the solutions of Eqs. (pT|) and ([22|) we obtained the values of Mg{T) and Dg{T) reported 
in Fig. [3l 
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FIG. 3: Left panel: Effective mass of the gluonic degrees of freedom in the quenched approximation as a 
function of the temperature. Right panel: Effective number of degrees of freedom of the gluonic quasiparticles 
in the quenched approximation as a function of the temperature. 



Both the effective number of degrees of freedom and the effective mass monotonically decrease 
with increasing temperature. Therefore, in our quasiparticle approach to the deconfined gluon 
plasma, gluons are propagating quasiparticles with a mass and an effective number of modes that 
decrease with increasing temperature. The reduction of the number of the gluonic modes with 
temperature is not new and is physically due to the fact that the spectrum of the gluons depends 
on the temperature. For momenta k < gT not only the transverse components are present, but also 
the longitudinal ones. On the other hand, at high momenta k » gT the longitudinal modes are 
not relevant [s^l because the corresponding pole in the propagator becomes exponentially small. 
Since the equation of state is dominated by particles with momenta ~ T one can expect that 
at low temperatures the longitudinal modes are excited, whereas in the high temperature regime 
the contribution of the longitudinal poles becomes negligible. In principle we should treat the 
longitudinal modes employing a different dispersion law, however for the sake of simplicity we 
have assumed that all modes are degenerate in mass. This is a first possible ingredient of a more 
complete dynamical description of the system that may lead also to avoid the divergent behavior 
of the effective number of degrees at the transition temperature. A second possible motivation for 
the sudden increase of the number of degrees of freedom close to Tc is that other quasiparticles, 
like glueballs, may be present (analogously to the fermionic correlated states that survive above 
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Tc). As a matter of fact it has recently been pointed out in Ref. 0| that the mass of scalar and 
pseudoscalar glueballs in the deconfined phase is much less than in the confined phase. Therefore 
one can expect that these states significantly contribute to the various thermodynamical quantities. 

A third possibility is that the dispersion relation in Eq. (|28p is strongly modified in the infrared 
region as in Ref. 291 ]. 

In any case, even with our very simple model, we obtain that at high temperatures the number 
of the gluonic modes is roughly 15, not far from the expected result 16, that corresponds to the 
number of transverse modes. 




H 
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FIG. 4: Gluonic pressure (left panel) and energy density (right panel) per degree of freedom as a function 
of the temperature. 



In Fig. S] we show the plot of the pressure pg and energy density ig per degree of freedom. 
Note that the energy density does not suddenly saturates at T ~ 1.2Tc as observed in the lattice 
simulations. The reason is that we have subtracted from the energy density the contribution due 
to the gluon condensate. 

A quantity that is relevant for the dynamical properties of the system and gives a different 
measure of the deviation from the conformal behavior is the velocity of sound, defined as, 



dp 
Ye 



(25) 



For a conformal symmetric system the velocity of sound squared is equal to 1/3 and therefore 
deviation from this numerical value indicate a breaking of conformal symmetry. 
In our case we can define the velocity of sound of the gluon quasiparticles by, 



den 



(26) 



which is clearly different from the total velocity of sound in Eq. (I25p . The two results are shown 
in Fig. m the dashed line (blue online) corresponds to Eq. ()25p , whereas the full line (red online) 
corresponds to Eq. (f26l) . 

The two velocity of sound are quite similar except for T ~ Tc. Unfortunately this is the region 
where the errors in lattice data are larger and one can assume that our result have at least a 20% 
error bar (for recent more accurate lattice data regarding the velocity of sound in pure SU(3) see 
[s^ ] ) . Note that the velocity of sound drives the hydrodynamical evolution of the system produced 
in heavy-ion collisions and in particular of the elliptic flow. Therefore a more detailed knowledge 
of the sound velocity would be helpful in understanding the dynamical properties of the system. 
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FIG. 5: (color online) The velocity of sound squared of the gluon plasma as a function of the temperature. 
The full line (red online) corresponds to the velocity of sound for the gluon quasiparticles. The dashed line 
(blue online) corresponds to the velocity of sound for the total system. 



III. ANALYSIS OF THE QUARK GLUON PLASMA 

Let us now consider the case of the Quark-Gluon Plasma where gluons as well as quarks are 
dynamical degrees of freedom. As in the previous Section we first consider the result of lattice 
simulations and show that the gluon condensate is not sufficient to explain the trace anomaly of 
the system. In analogy with Eq.(jlj) let us define cT''{T) the ratio between the chromo-electric 
gluon condensates at zero and finite temperature for finite quark mass niq. The pure gauge results, 
c^^(T) corresponds to the limit rriq — > oo. 



It turns out that in full QCD with a quark mass mq = 0.1 T, Ce (T) < cf^{T) for any T [251. l33l|. 
In the following we will refer to cT''{T) at niq = 0.1 T as Cg'^(T) and we will assume that it is not 
strongly dependent on the actual value of the quark mass. 

Following the same approach of SecHU in Fig. [6] we compare the unquenched results for (e — 
3p)/T^ with the gluon condensate contribution, given by the Equation §1 with < >o=i 0.02 



GeV and Cg'^(T) obtained by interpolating the lattice QCD results of Ref. [25 1. 

In general, the trace anomaly has a fermionic contribution proportional to < >. However, 
above the chiral symmetry restoration temperature, that we consider coincident with the decon- 
finement temperature, one has < qq >= 0. 
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FIG. 6: (Color online) Interaction measure (black) and gluon condensate (red) as a function of T/Tc- 



In full QCD we find a result similar to the one obtained in the previous Section for the gluon 
plasma: the gluon condensate is not the only contribution to the trace anomaly. Therefore we 
rewrite Eq. ([2]) as follows: 
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< >T= e - 3P =< >o - < >T +Ae , (27) 

where AG represents the contribution to the trace anomaly from gluonic, fermionic as well as 
bosonic degrees of freedom. 



A. Quasiparticle model for the Quark-Gluon Plasma 

We shall assume that above the confining temperature the system consists of fermionic quasi- 
particles with the quantum numbers of u, d and s quarks, gluon quasiparticles and in-medium 
mesons (or correlated pairs). 

Concerning the fermionic sector, we assume that the number of quark (antiquark) degrees of 
freedom Dq [Dq) is independent of the temperature and Dq = Dq = 18. This is motivated by the 
fact that the baryon chemical potential is close to zero and charm quarks are too massive to play 
a role in the temperature range that we are studying. The general form of the dispersion law for 
in-medium fermions can be written as, 



ooqik, T) = iOqik, T) = ^fc2 + m{k , T f + ^Rik, T) , (28) 

and has been studied in EH, where m and the self-energy have been evaluated taking into 
account the interaction of the quasiparticles in the medium. 

Here we will assume that for the relevant momenta, of the order of the thermal momentum, we 
can neglect m/k (see e.g. [19]) and treat the term S/j as a chiral invariant and temperature depen- 
dent effective mass M{T). Quarks are considered degenerate in mass because M[T) is generated 
by the interaction of the quarks with the medium and for light quarks this is independent of the 
value of the bare mass. From this assumptions the dispersion relation for u, d and s quasiparticles 
can be written as, 

U:q{k,T) = UJq{k,T) = k + M{T). (29) 

The structure of the in-medium correlated states as a function of temperature is not easily to 
valuate. These string like states may describe qq states as well as more exotic states [s^. We will 
consider that there are Di,{T) bosonic degrees of freedom and we expect that Dh(T) increases close 
to the Tc and decreases with increasing T, because for asymptotic values of the temperature the 
system is made up of quarks and gluons and Dh{T) must vanish. 

In the following we will neglect, as a first approximation, the effect on the thermodynamics 
quantities of the width of the bosonic states. Therefore we employ the dispersion law 



u;h{k,T) = ^k^ + M,{Ty, (30) 

where Mf,{T) is the in medium mass of the bosons. To evaluate the dependence of the numerical 
results on this parameter we have changed the mesons masses in the range M(T)—3M(T) obtaining 
variation of less than 15% of our results. 

Concerning the gluonic sector we will employ the results obtained in Section III A[ Assuming 
that the effective number of degrees of freedom and mass of gluons is given by the result obtained in 
the quenched case is clearly a rough approximation, because it assumes that the effect of dynamical 
quarks on the gluon quasiparticles is much smaller that the effect of the gluon medium. Nonetheless 
we expect that in the unquenched case the effective number of gluon degrees of freedom has a 
behavior similar to the one obtained in the quenched case. Indeed we expect that Dg{T) decreases 
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with increasing temperature going from 24 close to to 16 at high temperatures, for the reasons 
explained in Section IIIBI Therefore we will assume that Dg (T) is given by the same expression 
obtained in the quenched case. 

Regarding the temperature dependence of the in-medium mass of the gluons we expect in both 
the quenched and unquenched cases a similar behavior. The reason is that at asymptotic high 
temperatures the system has to become scale invariant and therefore the gluon mass must decrease 
with increasing temperature. In the following we will assume that the gluon mass in the unquenched 
case, Mg(T), is proportional to the gluon mass in the quenched case, Mg{T), that is shown in the 
left panel of Fig[3l However, since the numerical value of the gluon mass in the unquenched case 
can be different from the numerical value in the quenched case, we write 

Mg{T) = CMg{T) , (31) 

where c is a coefficient which parameterizes the variation of the gluon mass in the medium when 
dynamical quarks are present. In the numerical study we will consider the values c = 0.5, 1, 2. The 
corresponding dispersion law for gluon quasiparticles turns out to be 

CjgiK T) = ^k-^ + MgiTY . (32) 

Within these approximations, the expressions for the pressure, p, and energy density, e, are 
given by 

P = Pf +Pb+Pg (33) 



e 



ef + eb + ig + e^n , (34) 



where the subscripts /, b and g refer to fermionic, bosonic (in-medium mesons) and gluonic degrees 
of freedom respectively and where 



P^iT) =TA/^log(l±e--/^)±S 



d^k 

with i = f,b; the sign -|- (— ) refers to fermions (bosons), whereas the gluon contributions are given 
by 

d^k 



Pg{T) = TDg{T) J ^^og{l-e'^^/^r' 



and the contribution of the gluon condensate to the energy density is given by, 

Cn = ^<G2>o[l-cr(r)]. (37) 

In order to evaluate M[T) and D{y{T) we perform a simultaneous fit of the lattice data of pressure 
and energy density of 3 flavors quark matter with bare quark masses m = 0.4 T of Refs. 10, [s^ 
as a function of the temperature employing Eqs. ([55]) and ([5^ . For each value of the temperature 
we consider the central value of pressure and energy density of the lattice data. We will discuss 
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FIG. 7: (color online) Quasiparticle chiral mass (left panel) and effective number of bosonic degrees of 
freedom (right panel) as a function of the temperature for T ~ 1. — A.Tc- The gluon effective mass employed 
in the evaluation of both the quantities is given by Eq. ((3T|) with the coefficient c that takes the values 0.5 
for the dashed line (green online), 1 for the full line (black online) and 2 for the dot-dashed line (blue online) 
respectively. 



the dependence of our results on the numerical values of the lattice data in the following. In 



21[ we evaluated with an analogous method M(T) and Di){T) in the range 1.2 — 2.Tc, keeping 



the mass of the gluons and their effective number fixed; here we extend such a study to a larger 
temperature range and considering temperature dependent gluon mass and effective number of 
degrees of freedom. 

In Fig. [7|we show the quasiparticle chiral mass (left panel) and the effective number of bosonic 
degrees of freedom (right panel) as a function of the temperature for T ~ 1. — 4.Tc. The various 
line correspond to different values of the parameter c in Eq. (jSlh that is the ratio between the 
gluon mass in the quenched and unquenched cases. In the three cases considered the chiral mass 
rapidly decreases with increasing temperature. This behavior suggests that the breaking of scale 
invariance is mainly due to the gluon mass and the gluon condensate. 

Also the effective number of bosonic degrees of freedom rapidly decrease with increasing tem- 
perature. At temperature larger than 1.5 Tc, we find that less than 20 bosonic modes are excited. 
The divergence of Df,{T) for temperature close to Tc is due to the same numerical problem that 
we have discussed in the previous Section when we have presented the results regarding the tem- 
perature dependence of Dg{T). However in this case the exponential growth of Di, begins at a 
value of the temperature ~ 1.2 Tc- We have checked, that the numerical behavior below ~ 1.2 
is strongly dependent on the numerical value of the gluon condensate and of the energy density 
and pressure. On the other hand for temperatures larger than ~ 1.2 Tc our results are basically 
almost independent of the value of the lattice data. Considering different values of energy density, 
pressure or gluon condensate within the statistical error bars determines a variation of less than 
10% of our results. 

In Fig. [8] we show the contribution to the pressure and energy density due to the fermions (full 
black line), gluons (dashed red line) and correlated bosons (dot-dashed blue line) quasiparticles 
employing a gluonic mass given by Eq. (jSip with different values of c. 

In all the considered cases we find that at high temperature the dominant contribution to energy 
density and pressure is due to the fermionic modes. However at sufficiently low temperature 
the contributions of gluons and correlated pairs become relevant. In particular, the correlated 
states give a non negligible contribution to energy density and pressure for temperature ^ 1.5 Tc. 
Unfortunately, as we have already stressed, for temperatures smaller than 1.2 Tc our results for the 
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FIG. 8: Contributions to the pressure (left panels) and energy density (right panels) of fermionic (full black 
line), correlated states (dot-dashed blue line) and gluonic (dashed red line) states. The in- medium gluonic 
mass is given by Eq. (j31[) where we have employed different values of c. From top to bottom we have taken 
c = 0.5, c = 1 and c = 2 respectively. 

boson modes are not reliable. However, it is reasonable that the qualitative behavior that we find, 
with a peak in energy density and in the pressure at a certain value of the temperature, does not 
depend on our approximations and is not an artifact of the numerical errors. The reason is that 
close to Tc energy density and pressure must be zero, because of the normalization of the lattice 
data. On the other hand at large temperatures we expect that the system consists of quarks and 
gluons, the correlated pairs being melted and therefore cannot contribute to the thermodynamical 
quantities. Then Pb{T) and e;,(r) must have a maximum for a certain value of the temperature. 

In Fig. [9] we show the squared sound speed of the various quasiparticles for three different 
values of the parameter c defined in Eq. (I3ip . At high temperatures the fermionic quasiparticle 
have the largest velocity of sound, however at moderate temperature it is not clear which component 
is dominant. Note also that fermionic and bosonic quasiparticles seem to be "approximately 
conformal" for any value of c, in the sense that the corresponding velocity of sound approaches 
the conformal values already at temperatures of the order of 2Tc. In the case c = 2 and c = 1, 
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gluons seem to be more scale dependent. However when we reduce the gluon mass taking c = 0.5 
also gluons seem to approach the conformal limit at a temperature close to 2Tc. Regarding the 
correlated pairs, we have shown their behavior up to the temperature where Di, is different from 
zero. 




T/T^ T/T_, T/T^ 



FIG. 9: (color online) Sound velocity of the various components. Full line (black online) correspond to 
fermions, dashed line (red online) corresponds to gluons and dot-dashed line (blue online) corresponds to 
correlated pairs. From left to right the three panels correspond to the values c = 2,1,0.5 respectively for 
the in-medium gluonic mass given by Eq. (|3ip . 



IV. CONCLUSIONS AND OUTLOOKS 

In the present paper we have described the lattice QCD data of pressure and energy density in 
terms of a quasiparticle model. Such description of IQCD data is not new. However our results 
show that there are two dynamical ingredients, not taken into account in previous analyses, which 
play an important role: the difference between trace anomaly and gluon condensate and the survival 
of correlated, string like, pairs for T > Tc, both suggested by different lattice simulations. 

The pure SU(3) gluon plasma lattice results can be described in terms of gluon condensate 
and gluonic quasiparticles whose mass decreases with increasing temperature and whose effective 
number of degrees of freedom tends to the expected value, 16, for large T. 

Analogously, numerical simulations results for QCD with finite quark masses can be fitted by 
the gluon condensate plus a mixture of gluonic and fermionic quasiparticles and bosonic correlated 
pairs. 

Also in this case the effective number of bosonic degrees of freedom and the in-medium masses 
decrease with increasing temperature. At T ~ 1.5 Tc only the correlated pairs corresponding to 
the mesonic nonet survive and they completely disappear at T ~ 2rc. 

The temperature dependence of the sound velocity has been studied to give indications for the 
hydrodynamical models of relativistic heavy ion collisions. 

The present analysis can give partial answers to important problems. For example, lattice data 
show a deviations of the pressure and of the energy density of the system from the SB values of 
an ideal gas of quarks and gluons even at temperatures T > 2>Tc. At T = 3 Tc ~ 500 MeV the 
running coupling constant is large and the deviations from the SB behavior can only be partially 



understood with improved perturbative methods [371. l38l|. Our results show that at large T there 
are still effects due to quasiparticle masses that determine variation from the SB values, whereas the 
contributions of gluon condensate and of the correlated pairs turns out to be small. The dominant 
contribution to the deviation from the ideal gas results seems to be due to the mass of the gluons. 

Our conclusions follow by considering a temperature independent chromo-magnetic gluon con- 
densate that has been evaluated in lattice simulations only up to T ~ 2 Tc. A gluon condensate 
growing as T^ for large T suggested in Ref. (39|. would modify the whole picture of the system 
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for T » Tc- However it is not clear at which temperature the results of Ref. [39|] would be rele- 
vant. For this reason we have assumed a constant value of the chromo-magnetic gluon condensate, 
extrapolating lattice data up to 4Tc. 

Another limitation of the present study is the zero width approximation of the correlated pair. 
The effect of the width would be to modify the expression of energy density and pressure of 
the correlated pairs, but also to change in a self-consistent way the dispersion law of fermions 
1^ . 1^ . 4o| . These aspects and a more refined treatment of the gluon condensate will be analyzed 



in a forthcoming paper. 
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APPENDIX A: THERMODYNAMICS CONSISTENCY 

In order to describe the pressure and energy density of the Gluon Plasma in Section III Al and 
of the Quark-Gluon Plasma in Section IIII Al we have employed a quasiparticle expression for the 
various thermodynamical quantities. Since we have assumed that some parameters are temperature 
dependent the thermodynamics consistency of our results must be carefully checked. In general we 
can write the relation between pressure and energy density as follows, 

T^=P + e + Cr, (Al) 

where the correction Cr depends on the temperature. The thermodynamics consistency requires 
that 

Cr = 0. (A2) 

The correction to the thermodynamical relation for the Gluon Plasma is due to the temperature 
dependence of the parameters Mg{T) and Dg{T) and to the introduction of the gluon condensate 
and is given by, 

Cr{T) = Tm/-^^^ f dkk'- + W^P^ - '"on ■ (A3) 

dT ZTT^ J 1 — exp [ujg/T) tOg Dg dT 

Requiring that C^"(r) = 0, is equivalent to requiring that Dg satisfies the differential equation 

A'^ + BDg + C = ^, (A4) 

where the coefficients can be obtained from Eq. (|A3|1 and read 

A = Tpg (A5) 

B = tM„^27t' I dke (A6) 

^ dT J 1 - exp (LOg/T) LOg ^ ' 

C = -el^^ . (A7) 

Substituting the values of Dg{T) and of Mg{T), obtained fitting the lattice data of energy and 
pressure, in the previous equations, represents a self-consistent check of our result. We find that 
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the differential equation ()A4p is satisfied with a good accuracy and the corresponding correction 
C^" is of the same order of the error in the lattice data. Therefore in our method the differential 
equation (|A4p turns out to be selfconsistently satisfied within the range of values of the statistical 
errors in lattice data. 

We can treat the correction to the thermodynamical equation for the Quark-Gluon Plasma in a 
similar way. In this case such a correction is due to the temperature dependence of the parameters 
M(T), Df,{T), Mg{T), Dg{T) and to the introduction of the gluon condensate. In this case we 
obtain, 

dM 2Da f „,o 1 



1 + exp {u:q/T) 



+r4M— — - — - / dkk ; — — ;;^r?'6 

dT 27r2 y 1 - exp {cob/T) uJb Db dT 

- dMa Dg f r, 1 IT dDa , , 

Also in this case the thermodynamical consistency can be cast as a differential equation for the 
effective number of degrees of freedom. Assuming that the dependency of Dg on the temperature 
is the same as in the quenched case we obtain a differential equation for Db of the form 

A'^+BDb + C = Q, (A9) 
where the coefficients A, B and C are given by 

A = j dkk^\og[l-e^v{-^b/T)] (AlO) 

T [ 1 AMdM , , , 

B = n / dkk^ 1 — 77^ ^ (All) 

J -l + eyL^{uJblT) UJb dT ^ ' 

C = -T^ fdke (A12) 

2^2 J 1 + exp (ug/T) dT ^ ' 

^^-dMgDg f „,n 1 1 T dD„ , 

+TMg / dkk^ — — + - ^rL ■ (A13) 

^ dT 27r2 i 1 - exp {Cbg/T) Cbg Dg dT ^"'^ ^ ^ 

Substituting the values of Dg(T), Mg{T), Db{T) and M{T) obtained fitting the lattice data of 
energy and pressure, we find that the differential equation ()A9P is satisfied with a good accuracy. 
Also in this case the correction C^^ is of the same order of the error in the lattice data. 
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